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Using a chiral random matrix theory we can now derive the low energy partition functions and 
Dirac eigenvalue correlations of QCD with different chemical potentials for the dynamical and 
valence quarks. The results can also be extended to complex (and purely imaginary) chemical 
potential. We also discuss possible applications such as fitting to low energy constants and under- 
standing the phase diagrams of the full and partially quenched theories. 
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1. Introduction 

Since the work of Banks and Casher [[lj] it has been known that the low eigenvalues of the 
QCD Dirac operator are related to the breaking of chiral symmetry. Later Leutwyler and Smilga 
[^] showed that the low energy partition function for QCD was universal and could be determined 
directly from the finite volume chiral Lagrangian. This led to a set of sum rules for inverse powers 
of the Dirac eigenvalues. Shuryak and Verbaarschot [[J then discovered that these sum rules could 
also be obtained from a Random Matrix Theory (RMT) based on chiral symmetry. Thus the low 
eigenvalue spectrum of the RMT is also universal and agrees with QCD. By now this agreement 
has been extensively checked by lattice simulations. 

Since the RMT results contain an explicit volume dependence one use for them has been 
to extract infinite volume results from finite volume lattice simulations. This was first done for 
the chiral condensate [Q] since that was the only available parameter in the RMT. Since then it 
was realized that the low energy partition function of quenched QCD with a chemical potential 
also depends on the pion decay constant as does the eigenvalues The determination of exact 
results for the low eigenvalue distributions of quenched QCD with a chemical potential [^] has then 
provided one way to extract both constants in quenched lattice simulations [0]. The eigenvalue 
distributions for unquenched QCD are also know [|8|], but due to the complex action they are not 
practical for current simulations. One alternative that has been suggested is to use an imaginary 
isospin chemical potential [^]. Then the action is real and also a Hermitian eigensolver can still be 
used. However for the unquenched simulations this required generating lattices with an imaginary 
isospin chemical potential which, although it does not cause any complications, it is not generally 
done. 

The ideal case is to consider partially quenched ensembles where the dynamical simulation 
is performed at one chemical potential (preferably zero) and the eigenvalues are extracted using a 
different (nonzero) chemical potential. We thus consider observables of the form 

1 r Nf 
(0) = - d[A] OW(no)) e" s 'W n<fct[0(M/) + ™/] (1.1) 

A J 7=1 

where S g is the gauge action and ^(ju) is the Dirac operator with quark chemical potential /x. If the 
observable only depends on the low energy (zero momentum) modes then the results are universal 
and can be evaluated using a chiral effective theory such as RMT. The observables we consider 
here are correlations of the Dirac eigenvalues so that 




(1.2) 



with Xu{X) the complex eigenvalues of X. Partially quenched eigenvalue correlations with an imag- 
inary isospin chemical potential have recently been published [[K]] and details of the calculation for 
a real chemical potential will appear elsewhere [[IT]]. Here we present some details on the RMT 
used along with a sample of the results. 



2 



Partially quenched QCD with a chemical potential 



James C. Osborn 



2. Dirac eigenvalues from partially quenched partition functions 

One way to calculate correlations of Dirac eigenvalues is from the partially quenched partition 
functions 



Z NfM ({m f \m b },{ix f \^}) = / rf[A]e-*W ^ l ™ \ /J . (2.1) 



The partially quenched eigenvalue density for one dynamical flavor, for example, is given by flT2[ ] 
1 



pi(z,^o|mi,^i 



.(2.2) 

m b =z,m b =z,H b =IA> 



The presence of the extra (conjugate) flavors with negative /J. is required when the eigenvalue 
spectrum of the Dirac operator is complex. The necessary partition functions can then in principle 
be obtained from zero momentum chiral Lagrangians, however these are only known for certain 
cases. The purely fermionic partition function is given by ^ 

Z N , = [ dUtet{Uyzxv(\zVTvM{U + U^}-\F 2 VTv[U,Q][U\Q] J (2.3) 
Ju{lf f ) \2 4 J 

where £ and F are the tree level condensate and pion decay constant and M and Q are diagonal 
matrices of my and \if respectively. The bosonic partition function Zo,2 is also known [Q]. So far 
no direct evaluation of a partially quenched chiral Lagrangian has been performed. Instead the 
partially quenched partition functions can now be evaluated from RMT 



3. RMT for QCD with a chemical potential 

Here we will consider a variant of the RMT introduced in ||8|]. The Dirac matrix in a chiral 
basis is written as (for zero mass) 

9 f =( ° i"f*+bf®\ (31) 
■* S \ia f &+b f ®T J ■ { ' 

Here <1> and are complex (A^ + v) xN matrices with v the topological charge. We start with a 
more general form of the Dirac matrix where af = a(T,fif) and bf = b (T, jUy) are functions of the 
temperature and chemical potentials that we will determine next. The QCD partition function with 
Nf quark flavors can now be modeled as 

Nf 

Z= (d<5>d<d exp(-«A^Tr[<I> + <I> + @ t @]) ]~]det(% + my) . (3.2) 

The constant a sets the average eigenvalue spacing and will also be determined next. We could 
absorb it into a and b, but instead will choose a normalization such that a(T,n = 0) = 1 and 
b(T,n =0)=0. 
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The coefficients can be determined by comparing to the chiral Lagrangian at nonzero \x. The 



mapping follows very closely to that for jit = given in Q13|]. First the fermion determinants are all 
written in terms of Grassmann variables as 



det(%+m/) = J dy{ 2 d\j/{ 2 exp 



Then the Gaussian integration over the matrices <1> and & can be performed. The result is 



2 exp 



a f a g 



aN 



(3.4) 



If we consider only a baryon chemical potential such that all flavors have the same chemical po- 
tential \if = Hb, the coefficient of the four fermion term becomes [a 2 (r,/ifi) — b 2 (T , Hb)] / ccN . At 
zero temperature we expect that the partition function does not depend on the baryon chemical po- 
tential. We can ensure this in the RMT by setting a 2 (0, pL ) — b 2 (0, pL ) = 1 . For nonzero temperature 
this will not be the case and one can leave the functions a and b independent. As was mentioned 
previously JT4|], this could possibly be used to help map out the QCD phase diagram. 

A nonlinear sigma model is obtained by using a Hubbard-Stratonovitch transformation to 
break apart the four fermion terms and then integrate out the Grassmann variables. This gives 



doe 



det(a+M) A ' +v det(a T +P + Mf 



(3.5) 



with a aNfxNf complex matrix and P = Ao^A - Bo^B - such that P = at ]Xf = 0. We also 
define A, B and M to be the diagonal matrices of a/, bf and m/ respectively. For small m/ and pLf 
the determinants can be expanded as 

det(a+M) /v+v det(a t +P + M) A ' exp (A^Tr [ln(ff t (r) +Ma _1 + (M + P)o~^ 1 ]) . (3.6) 

At T = we can also use the relation A 2 — B 2 = 1 to expand A as A as 1 + B 2 /2 which we will use 
below. For my = and }Xf = the large ,/V saddle point is given by o^o = 1 /a. Thus the Goldstone 
manifold is simply the unitary group JT3|]. The low energy partition function is then 



dU det{U) v exp {NTr^faM{U + U^} - NTr[U, B] [t/ 1 ", B] 



(3.7) 



By matching terms with (2.3) we find that 



a 



b(T = 0,ii) 
a(T =0,ju) 



2N 

1 + \l 2 F 2 



V_ 
4N 



(3.8) 



With these definitions we then have a RMT that maps directly onto the zero momentum chiral 
Lagrangian. We note that the 0(pL 2 ) term in a above only contributes to the overall normalization 
of the partition function. We can therefore neglect it for most observables. However quantities like 
the partition function, particle number and related susceptibilities will depend on it. 
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Figure 1: Real part of the one flavor eigenvalue density for dynamical mass m = at different chemical 
potentials fi. The peaks have been clipped for better illustration at jj,F^/V = 2,4. 



4. Partially quenched eigenvalue correlations 

A nice feature of the above RMT is that the partition function can be rewritten directly in 



terms of the eigenvalues of the Dirac matrix ( p.l[ ). This was originally done for the unquenched 
model where the dynamical and valence chemical potentials are all the same It is also possible 
to do this for the partially quenched case, although the results become much more complicated. 
Due to space constraints we will only show plots of some of the results and save the details of the 



calculation for a future publication [11]. 

As a review we first show results for the one flavor eigenvalue density [|L O]. In figure |I] 
we plot the real part of the density for dynamical mass m = for different values of the common 
chemical potential }X. At /i = the eigenvalues are purely imaginary. For small pL they begin to 
spread out into the complex plane, but still show the characteristic oscillations along the imaginary 
axis due to chiral symmetry. For larger /I those oscillations disappear and a new set of oscilla- 
tions appear around the real axis. These oscillations grow exponentially with the volume and also 
change sign unlike the ones for small fi. They are not seen in the quenched density and are in turn 



responsible for producing a nonzero chiral condensate [ ]16[ ] . 

We now start with the one flavor density at \iF\/W = 2 and lower the dynamical chemical 
potential pL\ as the valence chemical is held fixed at }Xq = }X. This is show in figure ||. The peak 
that is seen near the real axis (which is clipped at fiiFy/V = 1.9) stays at about the same height 
as Hi varies. However the background density increases filling in the valley along the real axis 
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mFV 1/2 = 1-9 H 1 FV ,/2 = 1.5 




Figure 2: Real part of the partially quenched eigenvalue density for valence chemical potential ^LqF\/V = 2 
with dynamical mass m\ = and different dynamical chemical potentials ji\ . The peak has been clipped for 
better illustration at [l\F\/v = 1.9. 

until the oscillations are completely covered at }i\ = 0. At this point the density looks more like 
the quenched density shown in figure || than the one flavor result. The main differences between 
the quenched and partially quenched in this case are the overall scale and the relative size of the 
"bump" along the real axis. 

We can also extend these results to a complex chemical potential including a purely imaginary 
one. This could be used to test analytic continuation of quantities from imaginary to real chemical 
potentials. An imaginary chemical potential can also be used for fitting the low energy constants 
£ and F. However in this case it is better to consider the mixed eigenvalue correlation between a 
valence and dynamical eigenvalue since it is much more sensitive to the chemical potential. 

5. Conclusions 

We have shown results for the eigenvalue correlations of partially quenched QCD with differ- 
ent dynamical and valence chemical potentials. This could be useful for fitting low energy constants 
on existing lattices generated at zero chemical potential. It is also possible to generate the low en- 
ergy partition functions with different chemical potentials for each flavor from these results. This 
could then be applied, for example, to study three flavor QCD. Additionally the extension of the 
RMT to include coefficients which are general functions of temperature and chemical potential 
may be useful for mapping out the QCD phase diagram. 
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m fv" 2 .o 




Figure 3: Partially quenched (left) and quenched (right) eigenvalue density for valence chemical potential 
PLqF\/V — 2. The dynamical chemical potential for the partially quenched is pi\ = with dynamical mass 
m i = 0. Note the different scales. 
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